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SUMMARY 


The method of numerical integration for calculation of "beam, deflec- 
tion is presented in matrix form, to give it the advantages which are 
inherent in an influence -coefficient method. Both the scalar method of 
numerical integration and the influence -coefficient method may he improved 
by introducing weighting matrices. Only distributed loading is considered. 
Examples are presented to show that the use of weighting matrices reduces 
the calculation time required to obtain a desired degree of accuracy. 


INTRODUCTION 


Two methods for calculating small deflections of non-uniform beams 
are in general use: a method of scalar numerical integration and a method 

employing Influence coefficients. Either method may be expressed in matrix 
form to obtain maximum efficiency in the use of calculating machines. 

The method of scalar numerical integration for obtaining deflections 
gives a l l quantities that may be needed in stress analysis^ that is, shears, 
bending moments, section torques, curvatures, twists, slopes, and deflec- 
tions. (See, for examples, references 1 and 2.) In problems such as 
vibration and aeroelasticity where only the deflections are of direct 
importance, matrix methods are preferable since they provide a direct 
linear relation between loads and deflections and avoid the intermediate 
calculations required in the scalar method. The method of influence 
coefficients gives deflections due to concentrated loads through direct 
linear relations but requires first the calculation of the influence 
coefficients and is not suited without modifications to accurate deter- 
mination of deflections due to distributed loadings. 

A matrix method of numerical integration incorporating weighting 
matrices and especially formulated for obtaining deflections due to 
distributed bending or torsional loadings is a feature presented herein. 

This method is based upon an equivalence between distributed and 
concentrated loading obtained by assuming the loading curve to be a 



2 


MCA TN No. 1827 


series of parabolic arcs. The weighting matrices for "bending contain 
weighting numbers introduced in reference 1 and applied in reference 2. 

The weighting matrices for torsion have not been previously published. 
Weighting matrices permit the use of lower-order matrices than would be 
needed in c ommonly used procedures to obtain a desired accuracy. 

As an introduction to the weighting methods, the commonly used 
methods of numerical integration and influence coefficients for distributed 
loadings, based on step loading distributions, are formulated in matrix 
notation. The weighting methods described include, in addition to 
numerical integration, procedures for using weighting matrices in con- 
junction with influence coefficients. 

In all, four matrix methods of different degrees of accuracy are 
considered, and each method gives deflection as a linear function of 
distributed loading through an array of coefficients which might be called 
influence coefficients for distributed loading" or "transformation coeffi- 
cients" to distinguish them clearly from standard influence coefficients. 

Eesults of the application of three of the methods to a uniform 
cantilever beam are compared with exact solutions. 


SYMBOLS 


a 

a 

E 

G 

I 

J 

K 

*ij 

L 

m 


segmental area beneath torsional load curve 

segmental area beneath twist curve (increment of rotation) 

Young's modulus of elasticity 

shear modulus of elasticity 

moment of inertia of cross-sectional area 

mass moment of inertia, per unit of length, about the 
elastic axis 

torsion constant 

influence function 

influence coefficient 

length of beam 

bending moment 


P 


concentrated lateral load 



MCA TN No. 1827 


3 


p distributed lateral load 

Q concentrated torsional load 

5 static moment of mass, per unit of length, about 

the elastic axis 

T section torque 

4 distributed torsional load 

v beam shear 

x spanvd.se coordinate 

y bending deflection (translation) 

a curvature 

a concentrated curvature (bend) 

p slope of bending deflection curve 

6 Simpson's numbers 

(f> torsional deflection (rotation) 

9 twist (rate of rotation per unit length) 

X length of a segment 

'o> circular frequency of natural vibration 

p mass per unit of length 

| alternate spanwise coordinate 


NUMERICAL HfTEGKATION 


If full advantage is to be taken of the new procedures to be presented, 
all relations between variables must be expressed in matrix form- In order 
to provide a simple illustration of the formulation of the matrix equations, 
the well-known process of numerical Integration will be considered in some 
detail. A single typical scalar equation is written to illustrate the 
linear relationship being considered- The complete set of linear relations 
are then written in .expanded matrix form- The expanded matrix equations 
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are -written for a "beam span divided into four equal segments. The form 
of the matrices will he sufficiently general, however, to indicate the 
correct manner of extension to a larger number of segments. The matrix 
equations are also written in contracted symbolic form and may then be 
considered applicable to a beam with aDy number of segments. 


Torsion 

A cantilever beam with distributed torsional load is considered first 
as shown in figure 1(a) . The loading curve is considered to be replaced, 
for purposes of numerical integration, by a step function as shown in 
figure 1(b) . The step function is, in turn, replaced by a set of concen- 
trated torques as shown in figure 1(c)- The value of the concentrated 
torque is obtained by multiplying the ordinate to the load diagram by the 
width of the step. As a scalar example is given by 


% =Xa 3 


(1) 


In matrix form the concentrated torques are given by 


V 


1 0 0 0 0 


*1 

Qe 


0 2 0 0 0 


12 

q 3 

_ X 
2 

0 0 2 0 0 


"3 



0 0 0 2 0 


*4 



0 0 0 0 1 


*5 


( 2 ) 


The validity of this matrix equation is immediately seen from an applica- 
tion of the rcw-by-column rule for matrix multiplication. (See reference 3’) 
In contracted form this matrix equation may be written as 

M = | H H (3) 

In this equation and [qj are column matrices, or column vectors 

and [a] is the square matrix of coefficients. 

The section torques may now be computed. The section torque distribu- 
tion is a step diagram as shown in figure 1(d). As a scalar example T^ is 
given by 
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T 3 = ^ ^ (4) 

The matrix equation for section torques hecc!mes 


T 1 


1 0 0 0 0 


V 

T 2 


110 0 0 


^2 

T 3 

= 

1110 0 


% 

T 4 


11110 


Qi, 

T 5 


11111 


% 


In this equation the matrix of coefficients may "be regarded as an inte- 
grating matrix or a summation matrix. Throughout the paper there will 
appear , in various equations, two types of integrating matrices. The 
two Integrating matrices differ from each other in having values of zero 
or unity on the principal diagonal. The two matrices are defined as 
follows: 


0 0 0 0 0 
1 0 0 0 0 



110 0 0 
1110 0 
11110 


( 6 ) 


1 0 0 0 0 


110 0 0 



1110 0 
11110 
Ill'll 


(7) 


Equation (5) may now he written as follows: 


M ■ N M 
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The section torques may he expressed in terms of the applied load hy 
substituting from equation ( 3 ) in equation (8)j thus. 


>] “ I M H H 


The next step in the numerical integration is to divide the section 
torques hy the values of GJ at the corresponding stations to obtain 
values of twist 9 (change of rotation per unit of length). (See 
figures 1(d), 1(e), 1(f), and 1(g).) A typical equation for the third 
segment is given hy 



( 10 ) 


In matrix form the twists are given hy 


0 1 


1 

0 

0 

0 

0 

1 H 


T 1 

0 2 


0 

0 

0 

-C 

H 

O 


t 2 

0 3 

= 

0 0 l/GJ 3 0 0 


T 3 



000 i/gj^ 0 


*4 



H 

0 

0 

0 


t 5 


( 11 ) 


In contracted form this equation becomes. 


[ e] = [i/gj] [t] 

= [gj]" 1 [t] (12) 

The inversion of a diagonal matrix is discussed in reference 3- The twists 
may he obtained in terms of the applied load hy substituting from equation (5 
in equation (12) thus, 


(13) 
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The twist distribution, as shown in figure l(g), is a step diagram- 
This diagram may he integrated to obtain rotations. In scalar form 0^ 
is given by 


4 , 


= 10 ^ + ^- 0 )^ + 10 ^ 


(14) 


In matrix form the rotations are 


given by 


A 


11111 


*1 



0 1111 


0 2 

h 

= \ 

0 0 111 


0 3 

4k 


0 0 0 1 1 


*4 



0 0 0 0 1 


0 5 



— 




In contracted form equation (15) becomes 


(15) 


M - x H 1 H 


( 16 ) 


The matrix is the transpose of jS-^j - The transpose of a matrix 

is obtained by rotating the matrix about its principal diagonal. (See 
reference 3*) Substituting from equation ( 13 ) in equation ( 16 ) gives the 
final formula for rotations as follows: 


[ <d = % f 1] H 1 M M [1] (it) 


The square coefficient matrices in equation ( 17 ) may be multiplied 
together to obtain a single matrix for convenience in deflection calcula- 
tions- The resulting matrix j^C^j is dependent only upon the structural 

properties of the beam and is defined as follows: 


>] - v hi ' H’ 1 [ z J M (i8) 

Substituting from equation (l8) in equation (17) gives 


M ■ M H 


( 19 ) 
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The elements of the matrix j^C^j are coefficients which express a linear 

relationship "between the rotations and the ordinates to the distributed 
load curve. Since the numerical method that has "been used for obtaining 
equivalent concentrated loads is not highly accurate, a large number of 
stations along the span must be used to obtain accurate deflections. 


Bending 

The calculation of bending deflections by numerical integration 
involves the same procedures as those used in torsion. The matrix equations 
are therefore written generally in contracted form. The diagrams of 
figures 2(a), 2(b), and 2(c) show the concentrated loads to be related 

to the distributed load ordinates p^ by the equation 

M-IMM 

Direct summation of the concentrated loads gives the shears according to 
the following equation: 

(81) 

Substituting from equation (20) in equation (2l) gives 

CiKMHM (22) 


The diagrams of figures 2(d) and 2(e) indicate that the moment -curve 
ordinates may be obtained by direct summation of the areas beneath the 
shear curve. The moments are thus related to the shears by the equation 


m l 


0 • 
0 
0 
0 

°i 


ti 

mg 


10000 


y 2 

m3 

= X 

11000 


t 3 

m^ 


11100 



*5 


11110 


75 


( 23 ) 
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This equation is written in expanded form, to Ill ustrate the use of the 
summation matrix jjcjj . In contracted form equation (23) becomes 

H = x [ e °] [ t ] m 

Substitution of values from equation (22) in equation (24) gives 

H - x l M M M H <«> 

The moment values must now be divided by the corresponding El values 
which are illustrated in figure 2(f). The resulting diagram defines the 
curvature of the beam. These curvatures may be considered, for purposes 
of numerical integration, to be the loading on a conjugate beam (reference 4) 
as shown in figure 2(g). In matrix form the curvatures are given by 


H - M H 
= H ' 1 H 


( 26 ) 


The matrix [jElj 1 is a diagonal matrix. Substituting from equation (25) 
in equation (26) gives 






(27) 


The loading on the conjugate beam may now be converted to equivalent 
concentrated loads in the same manner used for the original beam loading} 
thus. 



( 28 ) 


The concentrated loads cc^, shown in figure 2(h), may also be regarded 
as concentrated curvatures- These concentrated curvatures may be visualized 
as bends in a broken line. Substituting from equation (27) in equation (28) 
gives 



( 29 ) 
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The slopes shown In figure 2(i), are computed by direct integration 

according to the following equation: 

M'H'M • (30> 

Substitution from, equation (29) in equation (30) gives 


[p] = tt [ s o] ' [ A J H " 1 M M M t ] < 3i) 

The deflections, as shown in figure 2(j), are obtained by an integra- 
tion of the slopes according to the following equation: 

W ■ x M ' W (32) 

Substitution of values from, equation (31) in equation (32) gives 

[r] =Tr[%.] [ E o] [ A ] [ EI ] 1 [ E J M MM (33) 


Equation (33) may be written in the following form for convenience in 
deflection calculations, as explained in the section on torsion: 


>] = [°b] H 


(3W 


where 


Just as in the torsion case, the matrix jc^J defines the relation between 
ordinates to a distributed load curve and the ordinates to a deflection 


curve. 


INFLUENCE COEFFICIENTS 


When influence coefficients are available, they may be used to determine 
deformations directly without the necessity for numerical integration. The 
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influence coefficients may "be computed "by the integration method of the 
section dealing with numerical integration or they may he obtained experi- 
mentally. The relationships between load and deflection are given by very 
simple matrix equations and are shown only in contracted form. 


Torsion 

Influence coefficients define a linear relationship between a finite 
number of concentrated loads and, a finite number of ordinates to the 
deflection curve. The matrix of torsional influence coefficients will be 

indicated by JjKL^j . This matrix is used to express the relation between 

rotations and concentrated torque loads as follows: 


M - [**] H 


The c omm on method of determining equivalent concentrated loads from a 
distributed load is to replace the distributed-load curve by a step curve 
as used previously. The equivalent of [qT] from equation ( 3 ) may therefore 
be substituted in equation ( 36 ) to obtain 


W'lWWC®] (37) 

If the influence coefficients are to be computed by the common method of 
numerical integration, a simple matrix formula may be used for this purpose. 
A comparison of equations (17) and (37) shows that the influence coeffi- 
cients are given by the following formula; 

[Kt] =x[Ei|' [gj]' 1 ^ (38) 

Equation (37) > when compared with equation (19) , provides a second 
definition of the matrix ^C-JJ , as follows: 


0] - | [„] [a] 


Bending 

The matrix of bending influence coefficients will be indicated by 

The relation between bending deflections and concentrated bending loads is 
given by the following equation: 
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>] - M [ p ] 


(bo) 


Th.e column vector of concentrated loads is replaced "by a vector of 
distritmted-load ordinates "by substituting from, equation (20) in eq.ua 
tion (4o)j thus. 


[7] - | p*] [*][>] 


(bl) 


A comparison of equations (33) and (4-1) indicates that the influence 
coefficients could he 'computed from the following formula: 

pb] - v [ s i] ' [ 2 o] ' M H _1 M H w 

Equation (4l) when compared with equation (34) provides a second definition 
of the matrix , as follows : 

L°>] - 2 N [ A ] (43) 

NUMERICAL INTEGRATION WITH WEIGHTING MATRICES 


The matrix methods of numerical integration described in the foregoing 
sections correspond to the commonly used scalar methods. In those methods, 
the manner of converting distributed loadings into concentrated loadings is 
rather arbitrary. Such conversions usually lead to appreciable errors in 
the deflections when only a smal 1 number of stations is used. 

The following sections show that the arbitrariness of the common 
methods can be largely removed by regarding the loading curves as series of 

parabolic arcs- With such an approach, the matrix Qa~| in the equations of 

the foregoing sections is replaced by matrices which are designated 
"weighting matrices. ' 


Torsion 


Consider the distributed loading curve of figure 3(a)* If* for 
example, the part of the curve between ordinates q,g and q.^ is assumed 
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to be a second- degree parabola defined by ordinates qr>, 
area between ordinates q.g and a is given by 


0 


and 


fy* the 


=2' “ h ( 5 «2 + % ' I4) 


(44) 


The derivation of this formula is shown in appendix A. Proceeding in like 
manner for the other areas permits each area to be written in terms of three 
ordinates, as follows: 


a l 


1 0 
0 

H 

00 


*1 

*2 


05 8-1 0 


<3.2 

a 3 

11 

00 5 8-1 


a 3 

a 4 


0 0-1 8 5 


14 

0 


0 0 0 0 0 


s 

_ __ 







m 


In contracted form, equation (45) becomes 


H • h 

and N is referred to as a weighting matrix. The areas H correspond 

to increments in section torque as indicated in figure 3(b). The section 
torques T at the five designated stations are related to the Increments 
in section torque by the equation 




(46) 


H= M H 


(47) 


Substitution of values from equation (46) in equation (47) gives 




(48) 


In the section on torsion dealing with numerical integration, a 
diagonal matrix of values of l/GJ i was introduced as shown in equation (11) 

The individual values of GJ^ were assumed to be the average values within 

a bay. This assumption corresponds to a replacement of the GJ diagram 
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by a step diagram as shown in figure 1(f). It is now necessary to make 
a slight modification of the definition of the diagonal elements of the 

matrix jjjjj . The elements of this matrix must now he defined as the 

values of GJj_ at the stations where the rotations are to he determined- 

These values are illustrated in figure 3(c)* With this revised definition 


of the matrix 
equation 



the distribution of twist (fig- 3(&)) is given hy the 


M-[«ra 


(49) 


or, upon substitution of values from equation (48) in equation (49), hy 


(50> 

In figure 3(4) y the areas a correspond to increments of rotation as 
indicated in figure 3(e)- These areas are computed on the assumption that 
the twist curve can he represented hy a series of parabolic arcs- For 
example, the area eig is found by the formula 


k ” k (" S l + + *3) 


(51) 


Equation (51) is similar to equation (44), hut the order of the coefficients 
is reversed. ThiB reversal is not necessary hut is used because it provides 
a convenience in the matrix equations which will he derived and discussed 
subsequently. The expanded form of the matrix equation relating increments 
in rotation to twists Is 


a l 


5 8-100 


V 

®2 


-l 8 500 


0 2 

a. 

x. 

0 

1 

H 

03 

oi 

o 



3 

12 



3 

% 


o 

0 

1 

H 

CO 

Ui 


e 4 

0 


0 0 0 0 0 


_ 9 5_ 




— 




( 52 ) 
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The elements in the matrices of equation (52) can he rearranged in such a 
manner as to provide a convenience in the matrix equations that are to he 
derived- All elements of the first two matrices must he moved downward one 
position to give the equation in the following form: 


0 


0 0 0 0 0 


0 1 

5 1 


5 8-100 


e 2 

Sg 

*121 

ii 

-1 8 500 


d 3 

*3 


o-i 850 


9 k 

— — 


l 

0 

o 

1 

I- 1 
00 
1 V. 


0 5_ 


(53) 


Although the elements of the column matrix on the left-hand side of the 
equation are out of their natural order, application of the row-hy- column 
rule for multiplication reveals that all of the linear relationships defined 
hy equation (52) are preserved in equation (53)- In contracted form, 
equation (53) is written as 


[ 5 ] = M w i]”H 


m 


in which 


[%] 


is a weighting matrix* The double prime is used to 


indicate a double transposition, or double rotation, of the matrix* The 
matrix is first rotated about its principal diagonal in the usual manner 
of transposition. The transposed matrix is then given a rotation about 
the secondary diagonal. The possibility of relating the vector |~aj to the 

vector by means of the matrix jj-fpj , which has already been introduced, 

arises because of the precautions that have been previously taken in making 
particular arrangements of the elements within the various matrices. 


Substituting from equation (50) in equation (5*0 gives 

[0 - m [<]" H’ 1 [ E °] H [0 


(55) 


The rotations 0 are obtained by a summation of the increments of 
rotation according to the formula 




( 56 ) 
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From equation (55) and. equation ( 56 ) the final formula for rotations is 
obtained, as follows: 



M ' W " H 


-1 




(57) 


A comparison of equation (57) with equation (19) indicates that the matrix 
of coefficients fc^l may he computed by the following formula: 



(58) 


Since this formula for C-fcl incorporates the weighting matrices, this 


def ini tion is expected to yield more accurate solutions than the definition 
given by equation (l8), which is based upon numerical integration without 


weighting numbers. In equation ( 58 ) the product 
seated by a single matrix [«} 




may be repre- 


H-[ E o][«l] < 59) 

! it r -1 

The product pjj can be obtained by a double transposition of mJ: 

[“]" = [ E o] ' [ w l] <&> 

The validity of equation (60) is proven in appen d ix B. Substitution of values 
from equations (59) and ( 60 ) in equation ( 58 ) giveB the following formula 


for 


[ C t : 




(61) 


n 

The matrices jjblJ and are standard universal matrices which are not 

dependent upon the properties of the beam or its loading- These matrices 
have therefore been computed and are given in appendix C for systems of 5> 
7, 9 , and 11 stations. The form of the linear relationship between the 
rotations and the ordinates to the distributed load curve remains as given 
by equation (19) • 
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Bending 

The conversion of a distributed "bending load of the type shown in 
figure 4(a) into a set of equivalent concentrated loads of the kind 
illustrated in figure 4(c) may "be accomplished in a manner analogous to 
the conversion of a distributed torsional load into increments of torque- 
The conversion to be used has been presented and illustrated in scalar 
form in reference 2. As in the section on torsion, the distributed 
loading curve is regarded as a series of second-degree parabolic arcs- 
The principle used in the conversion is indicated in figure 4(b) , which 
shows the load applied to a set of simply supported sub-beams that react 
on the cantilever beam at the five designated stations. The reactions of 
the sub-beams on the cantilever beam are the concentrated loads P- The 
statical equivalence of the concentrated loading and the distributed 
loading is restricted to bending moments at the five designated stations, 
but the effects at these stations are the only effects of direct interest. 

By appropriate integration, the following typical formulas for the equivalent 
concentrated loads are obtained: 

At end station 1 

P 1 = k ( Tp l + 6p 2 _ P 3 ) ^ 62) 

and at the intermediate station 2 

p 2 = £2 (^1 + 10p 2 + P 3 ) (63) 

In expanded matrix form, the complete expression of the equivalent 
concentrated loads in terms of the ordinates to the distributed loading 
curve is ! 



Written in contracted form, equation (64) becomes 

H = k H H (65) 
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and JW 2 J is referred to as a weighting matrix. 

The average shears in. the hays of the team, shown graphically in 
figure 4(d), are found by a summation of the concentrated loads according 
to the formula 




( 66 ) 


Substitution of values from, equation ( 65 ) in equation ( 66) gives 

[tsHHH <«> 

The bending moments at the five stations (fig. Me)) are obtained by inte- 
gration of the shears with the formula 


H - > CM 0] 


( 68 ) 


and, upon substitution of values from equation ( 67 ), equation (68) becomes 

H-sNKlMW <®> 

At this point, the reader is reminded that the bending-mcment diagram for 
the concentrated loads is a broken line whereas the bending-moment diagram 
for the distributed loading is a smooth curve- Both diagrams, however, have 
exactly the same set of ordinates at the five designated stations if the 
distributed loading curve is truly a second-degree curve- 

The distribution of curvatures a, shown in figure 4(g) , is obtained 
by dividing the true bending moments by the appropriate values of El ht 
all stations along the beam- This distribution gives a smooth curve/ to 
which the ordinates at the five designated stations are given by the equation 



( 70 ) 



from equation ( 69 ) in equation (70) gives 
2 


Substitution of 
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According to the conjugate team, theory, the diagram of curvatures is 
considered to te a loading on the conjugate team. By the use of formulas 
typified ty equations (62) and ( 63 ), the distributed loading on the conjugate 
team is then converted into equivalent concentrated loads 5. on the con- 
jugate team as shown in figure 4(h) . In contracted matrix form, these 
equi valent concentrated loads are given as 


DQ = k M H (t2> 

Substituting from equation (71) in equation (72) gives 

(T3) 

The average shears in the four hays of the conjugate team, diagrama- 
tically shown in figure 4(i), give the average slopes of the actual team. 
They are obtained ty a summation of the concentrated loads on the conjugate 
team with the f ormula 


M-H'p] 


Substitution from equation (73) in equation ( 7 I 4 .) gives 

The deflections of the actual team are equal to the tending moments in the 
conjugate team. Ordinates to the deflection curve of the actual team, 
shown in figure 4(j), are found ty integration of the slopes of the actual 
team according to the formula 

H = x [ 2 oJ'W (76) 

The final formula for deflections is obtained ty substituting from equa- 
tion (75) in equation ( 76 ) as follows: 


1 t 


>] - 5-76 Po] [ 2 l] [ W s] [ EI ] _1 [ E o] H [ w 2 ] [s_ 
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A comparison, of equation (77) with equation (3^-) Indi cates that the matrix 


may he def ine d as follows: 

A 


[C b ] 

[ c »] ■ 5 fs [ e °] ' [ E i]' &Q H _1 W Qa] Kl (78) 

Since this formula for includes weighting matrices, it is expected 

to yield more accurate deflections than the formula given by equation (35) 
obtained by numerical integration without weighting matrices. 

The formula for may be simplified by introducing the following 

definitions: 




The validity of equation ( 8 o) is proven in appendix B. Substituting 
equations (79) and (8o) in equation (78) gives the following formula 


for [C b j: 


from 


M- 5 4h"[ ei T 1 [ h ] 


m 


The matrices 



!! 

and J^nJ are standard matrices which are independent 


of the properties of the beam or its loading. They have been computed and 
are given in appendix C for systems of 5, 7 , 9 > and 11 stations. The form 
of the linear relationship between the deflections and the ordinates to 
the distributed-load curve remains as given by equation (3^) • 


OTFLTJENCE COEFFICIENTS WITH WEIGHTING MATRICES 


The use of influence coefficients measured on the structure may be 
desirable in the final stages of design as a check on preliminary calcula- 
tions or in order to account for the effects of items such as large 
discontinuities in structure and the restraint of warping of cross-sections* 
Weighting numbers can be introduced that will increase the accuracy of 
deflection calculations for distributed loading when influence coefficients 
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are used. Weighting matrices involving the appropriate ■weighting numbers 
are introduced in the following sections. 


Torsion 

In the method of numerical integration with weighting matrices, no 
consideration of concentrated torsional loads was necessary. It is therefore 
necessary to develop a concept of an equivalent concentrated torsional load 
for use with influence coefficients. This concept is developed in detail 
in appendix A and the result is stated here- If the distributed torsional- 
loading curve of figure 3( a )j for example, is truly a second-degree curve and 
if GJ is constant over the length of the "beam, the equivalent concentrated 
torsional loads are given exactly "by the formula 


Ql 


Qg 


% 

11 

% 


S 



7 6-1 

2 20 2 

0 2 20 

0 0 2 

0 0-1 


0 0 


*1 

0 0 


*2 

2 0 


*3 

20 2 


% 

6 7 


h 


( 82 ) 


In contracted form, equation (82) "becames 


M - k N H 


(83) 


If the equivalent concentrated torsional loads j^Qj are used in conjunc- 
tion with Influence coefficients, the deflection obtained will be exact 
for a uniform beam. The rotations are obtained by substitution of values 
from equation (83) in equation (36) as follows: 


M'sMHW (a) 

Equation ( 84) when compared with equation (19) provides another definition 
of the matrix [Ct] , as follows: 
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This definition of j^C^J involve a influence coefficients and weighting 

n umb ers should prove to he more accurate than the formula given by 
equation (39) which involved influence coefficients without weighting 
numbers* 


Bending 

When tending deflections due to a distributed load are computed with 
the use of influence coefficients, the method of equivalent concentrated 
loads gives more accurate results than the common method involving the 
step-function appr oxima tion* The equivalent concentrated loads as given 
by equation (65), however, are appropriate for the purpose of computing 
bending moments but not for the purpose of computing deflections when used 
with the influence coefficients. Certain advantages may consequently be 
obtained if a new approach is made to the deflection problem in order to 
develop weighting matrices for use in conjunction with influence coefficients* 

When an analytical influence function K(x, I) is known for the 
deflection at station x due to a unit load at Btation £, the deflection 
due to a loading p( 5) may be expressed by means of a definite integral 
in the following form (see reference 5* P* 2 66): 

y(x) =J~* K(x,l)p(£)d| (86) 


Equation (86) must be replaced by a system of linear algebraic equations 
by letting x and | take on a finite set of values x^ and corres- 

ponding to equally spaced stations. The evaluation of the definite integral 
can be performed accurately by using Simpson's rule* The value of y at 
station x^ may be indicated by y^ and the value of p at station I 

J 

by pj* The value of K(x,i) for x = Xj_ and I = £ j may be indicated 

as an influence coefficient kpj* Simpson's numbers may be indicated 

by e . A single scalar equation of the system is written for deflection 
J 

at station 3 as an example 


- k U 6 l p l + k 


B2 £ 2 P 2 


+ k„e p + fc-iAp. + 


( 87 ) 


From this equation the influence- coefficients may be seen to be formed into 
a square matrix and the weighting numbers to be formed into a diagonal 
matrix in order to express the system of equations in matrix form. These 
matrices appear as follows: 
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k U k 12 
^2! k 22 



( 88 ) 



e ± 0 0 0 0 

0 0 0 0 



0 0 0 6 ^ 0 

0 0 0 0 e 5 


(89) 


Simpson's numbers are well known (see reference 5, p* 5) and. may be used 
in equation (89) to obtain 


1 0 0 0 0 



0^000 


0 0 2 


0 0 


0 0 0 4 0 
0 0 0 0 1 


(90) 


The system of linear algebraic equations which is equivalent to 
equation (86) may now be expressed in the following matrix form: 


[ y HM MW 


(91) 


Equation (91) when compared with equation (34) provides another definition 
of the matrix as follows: 




(92) 
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Since this formula for 



includes -weighting numbers it should provide 


more accurate deflections than the formula in equation (43) • 


APPLICATIONS TO NATURAL VIBRATION PROBLEMS 


Equations (19) and (34) may be 
equation, as follows: 


formally stated in a single matrix 


C t 0 


q. 


0 

_° °b 


p 


7 


(93) 


The submatrices 


and 


Cb 


in equation (93) are the 


relating distributed loadings to deflections and have been 
in four ways with varying degrees of- accuracy- 


Bquare matrices 
defined previously 


The loads may be applied statically or dynamically. Natural vibration 
is an example of dynamic loading, the dynamic loads being the inertia 
loads which accompany the accelerations of the system during vibration. 

The cantilever beam under consideration vibrates naturally in any of an 
infinite number of modes with each of which is associated a frequency of 
vibration. If the deflections are small and in the elastic range, the 
vibration in each mode is a simple harmonic motion in which all parts of 
the beam oscillate about the position of static equilibrium in phase with 
each other and with the same frequency. The dynamic or inertia loads for 
any mode may be expressed in terms of the deflections and mass properties 
of the beam and the frequency as follows: 


1 

2 

s 

1 

0 


= 




p 


Q. 1 
CQ 

1 


7 


(94) 


In equation (94) , cd is the frequency of vibration, £ pj is a diagonal sub- 
matrix of values of mass per unit length, J^sJ is a diagonal sub ma trix of 

static moments of mass per unit length about the elastic axis of the beam, 

Lt 

H 


and 


-*mj 


is a diagonal submatrix of moments of inertia of mass per unit 
length a&out the elastic axis of the beam. The elements of the matrix Q 

are sometimes called the coupling terms, and when jjsj is different from 
zero, each natural mode contains both torsion and bending deflections 
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Combining equation (9*0 and equation (93) gives the complete equation for 
natural vibrations, as follows: 


C t 0 


•I m s 


$ 

1 

7 

i — 

0 

0 

l°L 


S P 


1 

% 2 

y 


(95) 


Equation (95) may be consolidated for practical use by carrying out the 
multiplication of the two square matrices on the left side. The resulting 
equation may be stated in contracted form as follows: 




7 

1 

0 

Y 







y 

0> 2 

y 


(96) 


The matrix 



depends only on the geometric, mass. 


and elastic properties 


of the cantilever beam and therefore remains constant for all modes of 
vibration. The solution of the vibration problem consists in finding sets 


of numerical values of the vector 


0 

y 


and associated frequency cd which 


satisfy equation ( 96 ). Iteration is a convenient method for obta ining the 
solution. Attention is drawn to the fact that, although the actual 
cantilever beam has an infinite number of natural modes of vibration, 
equation ( 96 ) can be used to calculate only a finite number of modes 
approximating with varying degrees of precision the exact modes in the 
Initial range of the frequency spectrum. Equation ( 96 ) determines n - 1 
modes, n being the number of designated stations on the beam- The preci- 
sion with which equation ( 96 ) determines the exact modes is greatest for the 
first mode, becoming less as the order of the mode increases. | 


COMPARISON OF METHODS 


In order to provide a simple comparison of the relative merits of the 
various methods of analysis that have been described, a uniform cantilever 
beam has been analyzed for both static and vibrational deflections. Exact 
values have been obtained by solving the well-known differential equations 
which govern the deflections. In the static analyses the loading is assumed 
to be of triangular shape with a value of zero at the tip and a maximum value 
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at the root of the team- In the vibration analyses the mass axis is assumed 
to coincide with the elastic axis of the "beam so that torsional and bending 
vibrations may be considered separately. 

The four methods of deflection analysis that have been described are 
as follows: 

(1) Numerical integration 

(2) Influence coefficients 

(3) Numerical integration with weighting matrices 

(h) Influence coefficients with weighting matrices 

Since the first method is well known and is usually the least accurate of the 
four methods, it has been omitted from consideration in this section. 
Calculations of deflections have been made by the other three methods. The 
first two natural frequencies and the static tip deflection have been 
computed for 3> and 7 stations and are recorded in table 1. The exact 
values obtained by solving the differential equations are also shown in 
table 1. The torsional deflections corresponding to static loading and 
vibrational motion are governed by the following two differential equations, 
respectively: 


GJ 



(97) 


GJ 


d 2 0 

dx 2 


+ pr cu 


V = 0 


(98) 


In these equations p is the mass per unit length of the beam, r is the 
radius of gyration of the mass, and cd is a natural frequency. The bending 
deflections corresponding to static loading and vibrational motion are 
governed by the following two differential equations, respectively: 


El 



= P 


(99) 


d^y p 

El — f - pen y = 0 

dx 4 


( 100 ) 
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The second method of calculation involve s the use of influence 
coefficients. The influence coefficients for a uniform, beam can be easily 
computed from simple formulas. These formulas constitute the influence 
function, or Green's function, for the problem. The influence function 
is a solution of the homogeneous differential equation for static deflec- 
tions. It gives the deflections due to a unit concentrated load. Since 
the function is discontinuous at the point of application of the load, 
it must be defined separately for the regions on either side of the load. 
If the deflection is to be determined at station x for a unit load at 
station | (see fig. 5 ) , the influence function for torsion is defined 
by the following two formulas: 


Kt(x, £) 

ii 

(x> I ) 

(101) 

K t (x,l) 

X 

“ GJ 

(x < 1) 

(102) 


The influence function for bending is defined by the following two equations: 



= H 
6ei 

(3x - 5) 

(x>S) 

(103) 

K-jjCx,!) 

s 2 

6ei 

(31 - x) 

(x<6) 

(ioi0 


Using the exact values, as shown in table 1, permits computation of 
the percentage error resulting from the use of each of the calculation 
methods for each number of stations. Graphs of this percentage error are 
shown in figures 6 to 11. In most cases the absolute value of the per- 
centage error decreases with an increase in the number of stations. In 
all cases the use of weighting matrices with either numerical integration 
or influence coefficients brings about an appreciable reduction in the 
percentage error. 

From practical considerations of economy in calculation effort, these 
graphs show that, for an allowable percentage error, the use of weighting 
matrices permits the deflection analysis to be made with a smaller number 
of stations. In order to illustrate this point an examination of the graphs 
has been made to determine the number of stations required to obta in 
satisfactory accuracy defined to allow the following percentage errors: 
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Static tip deflection 1 

First vibrational frequency 1 

Second vibrational frequency 2 

The number of stations required to give satisfactory accuracy, according 
to tbe above definition, bas been read from the graphs and recorded in 
table 2. From an ins pection of this table it is seen that, in five out of 
six deflection analyses, tbe use of influence coefficients alone would 
require more than seven stations to obtain satisfactory accuracy • When 
weighting numbers are used with influence coefficients, satisfactory 
accuracy is obtained in five out of six of the deflection analyses with 
less t han seven stations* When weighting numbers are used with numerical 
integration, satisfactory accuracy is obta in ed in all six deflection 
analyses with less than seven stations. 

The example that has been used for illustration purposes is a uniform 
beam and the relative percentage errors illustrated in figures 6 to 11 cannot 
be considered as strictly of general applicability. This indicates clearly 
the need for future research in studies of nonuniform beams. Future research 
must also deal with the development of more accurate weighting numbers and 
practical methods of analysis with concentrated loads. 


CONCLUDING REMAKES 


The advantage of an influence coefficient method of deflection analysis 
is that it provides a direct linear relation between the loading and the 
deflection in explicit form. The same advantage may be obtained in a numeri- 
cal integration process, employing beam stiffness properties, if the analysis 
is expressed in matrix form- The linear relationships for distributed 
loading have been developed. 

The accuracy of both the numerical integration and influence coeffi- 
cient methods can be improved by the introduction of weighting matrices. 
Consequently, for a desired degree of accuracy, smaller matrices may be 
used. This procedure results in an appreciable saving in calculation 
time since the computing work varies as the square of the order of the 
matrices . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., December 7, 19^8 
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APPENDIX A 

DERIVATIONS BASED ON THE PARABOLA 
Formula for Areas 


In the method of numerical integration with weighting matrices, 
increments of section torque and Increments of rotation correspond, 
respectively, to increments of area under the curve of distributed 
torsional loading and increments of area under the curve of twists- In 
the calculations, these curves are appr oximat ed by a series of second- 
degree parabolic arcs defined by groups of three ordinates- The formula 
for increments of area under a second-degree parabola are derived as 
follows: In figure 12, the ordinates f-]_, fg, and f^, separated by 

the distance X, represent a typical group of three ordinates to a loading 
curve or a curve of twists. With the coordinate system shown in figure 12, 
the equation of the second-degree curve defined by the three ordinates 
may be written in the form of Lagrange's interpolation formula (reference 6) 
as follows: 


f 1 X - x 2X 


f . f - * 

X + f 2 X X 


I 3 x X - x 

s’ r x 


(Al) 


The area a-^ between ordinates a* 1 * 3 - f 2 is found - "by integrating the 

function f between the limits x = 0 and x = X, and the result is 


a i -J 7 f - 5 ( 5f i * *2 - f s) 


(A2) 


Formulas for other areas may be obtained by increasing the subscripts in 
equation (A2). For example, the formula for area ag in figure 12 would 
have the form 


&2 


■50 


5fo + 8f Q - f. 


0 


(A3) 


The area a^ ma 7 also be computed by using equation (Al) and integrating 
from X to 2X- The resulting formula is 


a 2 = ^ (-fi + 8f 2 + 5f 3 ^ 


(Ak) 
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A formula of the type given "by equations (A2) or (A3) must he used at the 
left end of a heam- The type of formula shown by equation (Ah) must he 
used at the right end of a beam. Either type of formula may he used for 
intermediate segments of the span. 


Formulas for Equivalent Concentrated Torques 

The concept of equivalent concentrated torques, necessary for the 
use of weighting matrices with influence coefficients, is based on the 
following condition of equivalence: In the beam shown in figure 13, the 

set of rotations, at the five designated stations, due to the concentrated 
loading must be identical with the set of rotations due to the distributed 
loading. If this condition is fulfilled, the concentrated torques Q of 
figure 13 (b) must produce a set of increments of rotation in the bays of 
length 1 equal to those produced by the distributed torsional loading of 
figure 13 (a) . 

Consider first the increment of rotation between stations 1 and 2 
produced by the distributed loading. The section torque at a distance x 
from the tip is given by 

P z 

T x = q dx (A5) 

U 0 


The twist at the distance x is given by 


9 


x ~ 



_i_ p x 


dx 


(A 6) 


The increment of rotation between stations 1 and 2 is then given by 


f$i2 - 


PI 


fU 


0 


ax = 


fix 


UO 


GJ. 


xU 


q dx dx 


(A7) 


0 


Consider next the increment of rotation between stations 1 and 2 pro- 
duced by the concentrated loading. The section torque between stations 1 
and 2 is equal to Qp. The twist at the distance x between stations 1 
and 2 is given by 



(A8) 
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The increment of rotation between stations 1 and 2 is 



0 X 


0 


dx = Qi 


nx 


dx 


0 ^x 


(A9) 


Since the two expressions of equations (A7) and (A9) for ^ are "to 
he equal, the following equation is obtained 


PX 


Qi 


dx 

GJx 


OX 


fix 


uo 


GJ, 


q dx dx 


(A10) 


Equation (A10) may be integrated if the values of q and GJ Z are given 

as functions of x. The present consideration will be limited to the case 
of constant GJ; equation (A10) can therefore be simplified to the form 

Qq = ^ J 1 q dx dx (All) 

Equation (All) is found to apply also when GJ 1 b a step function. The 
variation of q is assumed to be given by the second-degree parabola 
defined by the ordinates q^, and q^- The expression for q is then 

obtained from equation (Al) by substituting q for 

!l 1 - x 21 - x , x 21 - x ^3 
9 2 1 1 92 1 1 2 


f, as follows: 


x 1 - x 
1 1 


(A12) 


With this expression for q, equation (All) becomes, after integration. 


Q 1 = 24 ( T<1 1 + ^ ^ (A13) 

The expression for Qg is. figure 13 (b) will now be derived- It is 
necessary to consider Qp as consisting of two parts; a part Qgq asso- 
ciated with the distributed loading between stations 1 and 2, and a 
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part Qgg associated with the distributed loading "between stations 2 and 3* 

The sum of and Q23 equals In figure 13 (h), if the sum of Q-j_ 

and equals the total loading between stations 1 and 2 in figure 13(a), 

the rotations at the five designated stations of the beam due to the loading 
between stations 1 and 2 in figure 13(a) will be equal to .the rotations 
due to the loads and Q2i* This consideration provides the definition 

of Q21, as follows: 


fl \ 


Q21 - 


u 


0 


q. dx - q-l 


(AUO 


The integral in equation (Alh) is the area under the q-curve between 
stations 1 'and 2 and in accordance with equation (A 2 ) is given by 


D X 
0 


4 - h ( 5l i + • 43) 


(A15) 


Substituting from, equations (A 13 ) and (A 15 ) in equation (Alh) gives 


^21 ~ 2 k ( 3<1 1 + 10q -2 " ^3) 


(A 16 ) 


The concentrated load Q23 bears the same relation to the distributed 

loading between stations 2 and 3 us the concentrated load Qj bears to 

the distributed loading between stations 1 and 2; this relation gives that 
part of the total increment of rotation between stations 2 and 3 which is 
due to the distributed loading between stations 2 and 3 * The defining 
equation for similar to equation (All) is, as follows: 



(A 17 ) 


Substitution of values from equation (A 12 ) in equation (A 17 ) gives, after 
integration. 


^23 = £k 1 + 10 ^2 + 31 ^) 


(A 18 ) 



5K 


MCA TN Ho. 1827 


33 


The sum of equations (Al6) and (Al8) gives 


02 


\ 


= 15 




(A19) 


Formulas for the equivalent concentrated loads at the other stations may 
he derived hy proceeding in the manner used to derive equations (A13) and 

(A19) • 
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APPENDIX B 

DO UBLE TRANSPOSITION OF MATRICES 


A matrix j^A^j is transposed "by rotating it about its principal 

diagonal. The tr ans posed matrix is indicated with a prime as A • 

1 — — 1 » L J 


If the matrix 




is now rotated about the secondary diagonal the 

r.i" 


new matrix may be indicated with a double prime as AJ . The matrix 

may be said to be doubly transposed* These transposition processes 
may be illustrated by the following three equations: 




“tl 

*12 

^3 

M- 

^1 

**22 

^3 


> 

a 32 

a 33_ 


"^li 

®21 

V 

[a] . 

a i 2 

®22 

a 32 


> 

^3 

& 33_ 


a 33 

a 32 

a 31 

[A]" - 

^3 

a 22 

®21 


> 

S 12 

a ll_ 


(Bl) 


(B2) 


(B3) 


The object of the present appendix is to prove the validity of the 


formulas for 


For 


^M~j" and j^Nj" as given i n the main body of the paper. 

this purpose the several types of symmetry which matrices may have must 
first be considered. If a matrix is equal to its transpose, the matrix 
is said to be symmetrical. This condition might be referred to more 
specifically as a principal symmetry. A matrix that remains unchanged 
after rotation about its secondary diagonal might be said to have secondary 
symmetry. If a matrix remains unchanged after a double transposition. 
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the matrix is described as being centrosymmetrical (reference 7)* These 
three types of symmetry may "be expressed by equations relating the matrices 

[4 W, and [A] as follows: 


Principal symmetry: 


Central Symmetry: 


OH a]' 

Secondary symmetry: 

[a]’ =[aJ" (B5) 

Central Symmetry: 

[ A ] = [a]" (B6) 

A relationship show ing that any matrix j^A~J" can ^> e obtained from the 
matrix by simple matrix multiplications is now necessary. A matrix 

must be introduced according to the following definition: 


?]■ 


0 0 1 


0 10 


10 0 


This matrix has values of un i ty in its secondary diagonal and zeros else- 
where. If a matrix J^aJ is premultiplied by the procedure merely 

interchanges the rows of j~Ajj thus. 


a 31 a 32 a 33 


_j][a]- ag-L 

a 22 ^3 


11 a 12 a 13 


If a matrix |~aJ is postmultiplied by QlJ, the procedure merely inter- 
changes the columns of [| aJ j thus, 


a 13 a 12 a ll 


Q oo &oo Q.* 


a 33 a 32 a 31 
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If a matrix 



Is 


the procedure “brings 
and vertical central 


premultiplied by 



and postmultiplied by 



about a double rotation of tbe matrix about horizontal 
axes. Tbe result is as follows: 


HMH 


*33 

"32 

"31 

*23 

a 22 

*21 

"l3 

*12 

*11 


(BIO) 


If equation (BIO) is compared with equation (B 3 ) , the expanded matrix 
on the right-hand- side is found to be exactly the same in both equations 
and therefore 


[ A ]" -DQMH 

An unusual but important property of the matrix 



(Bll) 


must now be noted. 


The matrix is its own reciprocal, which may be expressed by the following 
equation: 


[ j T=[ j ][ j ]=L i ] ^ 


In this equation 



is the identity matrix. 


Development of the desired formulas is now possible, 
case will be considered first. According to equation (59) 
is defined by the formula 


The 

the 


torsion 

matrix 



M ° [ z °] M (E13) 

According to equation (Bll) , the doubly -transposed matrix may be determined 
from the equation 


M" -M W M M 


(B14) 


The identity matrix may be inserted as a factor at any point on either side 
of an equation without changing the value of either side. Equation (Bl4) 
may therefore be written in the following form: 


H" - M M H H W 
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Equation (B12) may now "be combined with equation (B15) to obtain 


H" =[ j ] H HH H H 


(B16) 


Formula (Bll) applied in equation (Bl6) gives 



(B17) 


An interesting feature of the process of double transposition can be noted 
at this point- Comparing equations (B13) and (B17) indicates that the double 
transposition of a product of matrices can be obtained by a double transposi- 
tion of the individual matrices without changing their order. 


It must now be noted that 
diagonal. Consequently, 



is symmetrical about its secondary 



(B18) 


Substitution of values from equation (Bl8) in equation (B17) gives the 
following formula: 

= (KL9) 



This formula was 


The bending 


the matrix 


H 


previously given in the main text as equation ( 60) . 

case is now to be considered. According to equation (79) 
is defined by the formula 


>] - [*Q [\] [»2 


(B20) 


Employing the rule that has just been developed permits both sides of 
this equation to be doubly transposed to give 


[”]"*■’ jfo]" M '' N " (B21) 

The matrices jpoj and ' hoth have secondary symmetry so that 

equation (B2l) may be written as 
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>oT M’ [»3" 


(B22) 


A consideration of "the expanded form of the matrix jjfgj ’ as s ^ Lown ^- n 
equation (64) , indicates that this matrix has central symmetry } hence 



(B23) 


Substitution of values of equation (B 23 ) in equation (B22) gives the 
following formula: 



(B24) 


This formula was given in the main text as equation (8o). 
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APEEHDH C 


STAHMKD MATRICES 


Tiie matrices H and [mJ serve to weight and integrate the 

torsional loading curve and the twist curve, respectively. The 

matrices |^hJ and serve to weight and perform a double integra- 

tion of the load curve and the curvature curve in "bending, respectively. 
These matrices are standard matrices which may he tabulated and used 
for various cantilever-beam analyses. The matrices must be computed for 
each different number of subdivisions of the Bpan. The order of each 
matrix will be one greater than the number of segments into which the 
span is divided. Equation (59) gives the following formula for the 


matrix 


M : 


M ■ M M 

The formula for the matrix equation (79), is 

H ■ M H M 


The matrices and £nJ 


are obtained from 


H by 


double transposition. Equations (Cl) and (C2) have been used for 
computing the matrices f M~1 and r-i. 
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The matrices of fifth order are as follows: 


Fifth Order Matrices 


[m [m 


~0 

0 

0 

0 

O' 


~k 

15 

11 

13 

5 ~ 

5 

8 

-1 

0 

0 


-1 

7 

12 

13 

5 

5 

13 

7 

-1 

0 


0 

-1 

7 

13 

5 

5 

13 

12 

7 

-1 


0 

0 

-1 

8 

5 

_5 

13 

11 

15 

4 _ 


_0 

0 

0 

0 

o_ 


M M 


0 

0 

0 

0 

0~ 


~2 

21 

kk 

88 

35 " 

7 

6 

-1 

0 

0 


0 

2 

21 

60 

25 

16 

32 

0 

0 

0 


0 

0 

0 

32 

16 

25 

60 

21 

2 

0 


0 

0 

-1 

6 

7 

2 k 

88 

kh 

2 k 

2 


_0 

0 

0 

0 

0 


It is seen in the foregoing examples that the matrices 
are obtained from J^mJ and by a double transposition- 


and 

Since these 


examples illustrate the transposition process clearly, the higher order 
matrices will not be shown in transposed form- The higher order matrices 
are as follows: 
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Seventh. Order Matrices 




1 

M 




0 

0 

0 

0 

0 

0 

5" 

5 

8 

-1 

0 

0 

0 

0 

5 

13 

7 

-1 

0 

0 

0 

5 

13 

12 

7 

-1 

0 

0 

5 

13 

12 

12 

7 

-1 

0 

5 

13 

12 

12 

12 

7 

-1 

5 

13 

12 

12 

11 

15 

4 




>] 




~~0 

0 

0 

0 

0 

0 

0 

7 

6 

-1 

0 

0 

0 

0 

16 

32 

0 

0 

0 

0 

0 

25 

60 

21 

2 

0 

0 

0 

34 

88 

44 

24 

2 

0 

0 

43 

116 

67 

48 

24 

2 

0 

52 

144 

90 

72 

48 

24 

2 
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Ninth. Order Matrices 



4f 


0 

0 

0 

0 

0 

0 

0 

0 

0 

5 

8 

-1 

0 

0 

0 

0 

0 

0 

5 

13 

7 

-1 

0 

0 

0 

0 

0 

5 

13 

12 

7 

-1 

0 

0 

0 

0 

5 

13 

12 

12 

t 

-1 

0 

0 

0 

5 
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0 
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0 

0 
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0 

0 

0 

0 

0 

0 
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25 

60 

21 

2 

0 

0 

0 

0 

0 

34 

88 

44 

24 

2 

0 

0 

0 

0 
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67 
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24 

2 

0 

0 

0 

52 
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90 

72 

48 

24 

2 

0 

0 

61 
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113 

96 

72 

48 

24 

2 

0 

70 
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96 

72 

48 

24 
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Eleventh Order Matrices 
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TABLE I 

HATHRAL BKEQHEKCHS AKD STATIC TIP EEFLECTIOKS 



Torsion 

Bending 


First 

frequency, 

Second 

frequency. 

Static tip 
deflection, 

First 

frequency. 

Second 

frequency, 

Static tip 
deflection. 

Humber of 
stations 

■ 4 ? 


IL 2 
jn. 

GJ 

d M- 

1 T 2_ 

L ^r 


JpiA 

“2 VEI 






(a) 

■ 


mm 


V. 


Influence-coefficient method 



3 

1-530 

3-11 

0.125 

3-16 

13-70 

0-0261 

5 

1.560 

4.38 

.156 

3-if2 

20-5if 

-0316 

7 

1-567 

If. 57 

.162 

3-47 

21.29 

.0326 

Weighted-influence -coefficient method 

3 

1.575 

5-39 

0.1 er 

3-56 

15.63 

0.0347 

5 

1-571 

If- 73 

.161 

3.52 

22.8o 

•0334 

7 

1-571 

lf-72 

.167 

3-52 

22.08 

•0334 

Weighted-integration method 

3 

1.582 

7-59 

o.ier 

3-58 

Ilf. 58 

0.0347 

5 

1.573 

if. 

en. 

.167 

3-52 

22.08 

-0334 

7 

1.572 

if. 74 

.167 

3.52 

22.26 

•0334 

Exact values 


1-571 

if. 71 

0.l£7 

3-52 

22.03 

0.0333 




^oad distribution triangular with intensity q^ at root- 
^Load distribution triangular with intensity p^ at root- 
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TABLE II 

STATIONS REQUIRED FOR SATISFACTORY ACCURACY 


Quantity 

Influence 

coefficients 

Weighted. 

Influence 

coefficients 

Weighted. 

integration 

Torsion 

°^L 

5 

3 

3 

“2 

>7 

5 

5 


>7 

3 

3 

Bending 


>7 

4 

4 

“2 

>7 

7 

5 

y l 

>7 

5 

5 
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Cc) Equivalent 

concentrated loads. (g") Twist. 



Figure l.~ Torsional deflections obtained by 

step functions. 



(cl) Shear. 


(i) Slope. 



(e) Moment. 


y, 



(j) Deflection. 



Figure 2.- Bending deflections ob'tained by 

step functions. 









(c) Stiffness. 



(d) Twist. 


fe) Rotation. 



Figure 3.- Torsional deflections obtained by 

parabolic arcs. 
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(b) Load acting 
through sub- beams. 



(c) Equivalent 
concentrated loads. 



(d) Shear. 


m 

5 



(e) Moment. 
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(g)Curvature. 



(h) Equivalent load on 
conjugate beam. 



(1) Slope. 



ef le c+ions 


obtained by 


i c arc s . 
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(a 1 ) Beam with 
unit torque . 



(b) Rotation. 



(c) Beam wi+b 
unit force. 



Figure 5 .- Deflections of a cantilever beam 

wi+h uni t loads. 
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Number of s+a+ions 



Figure 6.- Comparison of errors in first 

torsional frequency. 
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Weighted 
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coefficients 


Influence 
coeff icients 
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Number of stations 


Figure 7. - Comparison of errors in second 

torsional frequency- 
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Figure 8. - Comparison of errors in stafic 

torsional tip deflection. 











56 


NACA TN No. 1827 



O 1 2. 3 4- 5 6 7 8 


Number of stations 



Figure 10. - Comparison of errors in second 

bending frequency. 
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Number of 


stations 



Figure i 1.- Comparison of errors in static 

bending tip deflection. 





MCA TN No. 1827 




Figure 12.- Areas under a parabolic curve. 



. . Q . £—Q- 

Q4. Q3 Glj Q.| 


(a) Torsional 
loading. 


(b) Equivalent 
concentrated 
torques. 


Figure 13.- Equivalent concen+ra+ed 
torsional loading. 
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